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ABSTRACT 

We present photometry of the GO star HAT-P-1 during six transits of its close- 
in giant planet, and we refine the estimates of the system parameters. Relative to 
Jupiter's properties, HAT-P-lb is 1.20 ±0.05 times larger and its surface gravity 
is 2.7 ±0.2 times weaker. Although it remains the case that HAT-P-lb is among 
the least dense of the known sample of transiting exoplanets, its properties are 
in accord with previously published models of strongly irradiated, coreless, solar- 
composition giant planets. The times of the transits have a typical accuracy of 
1 min and do not depart significantly from a constant period. 

Subject headings: planetary systems — stars: individual (HAT-P-1, ADS 16402B) 

1. Introduction 

More than 12 years have passed since the surprising discovery of "hot Jupiters": giant 
planets around Sun-like stars with orbital periods smaller than ~10 days (Mayor & Queloz 
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1995; Butler et al. 1997). These objects, of which about 50 are known (see, e.g., Butler et 
al. 2006), probably formed at larger orbital distances and migrated inwards through processes 
that are not yet fully understood. Hot Jupiters are also interesting because they are more 
likely to transit their parent stars than more distant planets. Transits are highly prized 
because they permit the determination of the planetary radius and mass (Henry et al. 2000, 
Charbonneau et al. 2000), the infrared planetary spectrum (Richardson et al. 2007, Grillmair 
et al. 2007) and longitudinal brightness profile (Knutson et al. 2007a), the stellar obliquity 
(Winn et al. 2007a), and much more. This helps to explain why so many groups around the 
world are conducting wide-field photometric surveys for planetary transits. Over a dozen 
cases of transiting exoplanets have been identified in this manner, with the rest having been 
found as a by-product of Doppler planet surveys (see Charbonneau et al. 2006 for a recent 
review) . 

Recently, the Hungarian-made Automated Telescope Network (HATNet) announced the 
discovery of HAT-P-lb, a giant planet that orbits one member of a G0/G0 stellar binary 
(Bakos et al. 2007). This planet was notable for being among the largest and least dense of 
all the planets for which such measurements have been made — both inside and outside the 
Solar system — and therefore an interesting test case for models of planetary atmospheres 
and interiors. 

However, while the data presented by Bakos et al. (2007) was certainly good enough to 
clinch the case for planethood and to provide useful estimates of the system parameters, it 
is possible and desirable to improve the accuracy of those parameters with repeated, high- 
precision, ground-based transit photometry. This is one goal of our Transit Light Curve 
(TLC) project, which has been described at greater length elsewhere (see, e.g., Holman et 
al. 2006, Winn et al. 2007b). 

This paper presents our results for the HAT-P-1 system. The next section describes the 
observations. In § 3, we describe the parameteric model that was fitted to the data, and 
in § 4 we present the results for the planetary, stellar, and orbital parameters, as well as 
the transit times. At the end of this paper we discuss the significance of the refined radius 
measurement. 



2. Observations and Data Reduction 

Our observations took place in late 2006, using telescopes at three different observatories. 
We observed 6 distinct transits and produced 7 independent light curves. 

We observed the transits of UT 2006 Sep 18, Sep 27, and Oct 6 with the 1.2m telescope 
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at the Fred L. Whipple Observatory on Mt. Hopkins, Arizona. We used the 4096 2 KeplerCam 
CCD, which has a 23(1 x 23f 1 field of view. We employed 2x2 binning, giving a scale of 
0'/68 per binned pixel, a readout and setup time of 11 s, and a typical readout noise of 7 e" 
per binned pixel. We used a Sloan z filter in order to minimize the effect of atmospheric 
extinction on the relative photometry, and the effects of stellar limb darkening on the transit 
light curve. We kept the image registration as constant as possible. We also obtained dome- 
flat and bias exposures at the beginning and the end of each night. On Sep 18, the sky 
conditions were photometric and the seeing varied from 1"7 to 2'.'1. We used an exposure 
time of 15 s. The night of Sep 27 began with patchy clouds and large transparency variations, 
but the rest of the night was clear. The seeing varied between 1"7 and 2'.'4, and we again 
used an exposure time of 15 s. Most of Oct 6 was lost to clouds, although we did manage to 
observe the egress in l'/8 seeing, using an exposure time of 10 s. 

We observed the transits of UT 2006 Sep 1, UT 2006 Sep 10 and UT 2006 Sep 18 
using the Nickel lm telescope at Lick Observatory on Mt. Hamilton, California. The only 
night when a complete transit could be observed was Sep 18. We used the Dewar ^2 direct 
imaging detector, which is a 2048 2 Lawrence Labs CCD with a 6(1 x 6(1 field of view. For our 
observations we used 2x2 binning (0'/36 per binned pixel), and read out only a 1450 x 500 
pixel subregion of the chip to decrease the readout time. Setup and readout time took about 
10 seconds per exposure, with a typical read noise of 11.8 e~ per binned pixel. We observed 
through a "Gunn Z n filter (Pinfield et al. 1997). To draw out the exposure time and to 
spread the light from stars over more pixels, we defocused the telescope until the stellar 
images had a full-width at half-maximum (FWHM) of about 6 pixels. The exposure time 
ranged from 10 to 40 seconds, depending on seeing and transparency. All nights were fairly 
clear with l'.'O-l'.^ seeing. On Sep 18, near the transit midpoint, the star passed within a few 
degrees of the zenith and autoguiding failed. The data from that time period were excised. 

We observed the transits of UT 2006 Sep 14, UT 2006 Nov 20, and UT 2006 Nov 29 
using the lm telescope at Wise Observatory, in Israel. We used a Tektronix 1024 2 back- 
illuminated CCD detector, giving a pixel scale of C/7 and a field of view of 11(9 x 1 1(9. We 
observed through a Johnson / filter, the reddest optical band available on the camera. On 
Sep 14 and Nov 20, weather conditions were poor, with patchy clouds. Because the data from 
those nights were of much lower quality than the other data presented in this paper, in what 
follows we describe only the data from 2006 Nov 29. The night was not photometric, and 
the measured stellar fluxes varied by about 20% over the course of the night. The exposure 
time was 15 s, and the FWHM of stellar images was about l'/8 (2.5 pixels). 
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We used standard IRAFlj procedures for the overscan correction, trimming, bias subtrac- 
tion, and flat-field division. We performed aperture photometry of HAT-P-1 and 4-8 nearby 
stars, depending on the telescope. The sum of the fluxes of the comparison stars was taken 
to be the comparison signal. The light curve of HAT-P-1 was divided by the comparison 
signal, and then divided by a constant to give a unit mean flux outside of transit. 

We then assessed residual systematic effects by examining the correlation between the 
out-of-transit flux and some external variables: time, airmass, the shape parameters of stellar 
images, and the pixel position of HAT-P-1. For the FLWO data, the flux variations were 
most strongly correlated with airmass; for the Lick data, the strongest correlations were 
with the pixel coordinates, especially the row number; and for the Wise data, there were 
correlations with both airmass and the FWHM of stellar images (which were themselves 
strongly correlated). We solved for the zero point and slope of the strongest correlation as 
part of the fitting process described in the next section. 

Figures 1 and 2 show the final light curves. The bottom panel of Fig. 2 is a phase- 
averaged composite of the 3 best light curves. Table 1 provides the final photometry, after 
correction of the residual systematic effects. 



3. Determination of System Parameters 

To estimate the planetary, stellar, and orbital parameters, and the times of transit, we 
fitted a parameterized model to the photometric data. The model and the fitting method 
were similar to those described in previous TLC papers (see, e.g., Holman et al. 2006, Winn 
et al. 2007b). It is based on a circular orbit@ of a star (with mass M+ and radius R+) and a 
planet (M p , R p ) about their center of mass, inclined by an angle i relative to the sky plane. 
Because one of our goals was to measure the individual transit times, we allowed each transit 
to have an independent value of T c , the transit midpoint, rather than forcing them to be 
separated by regular intervals. 



1 The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 

2 This is our default assumption in the absence of clear evidence for an eccentric orbit. Although the 
orbital fit of Bakos et al. (2007) yielded the formal result e = 0.09 ± 0.02, we regard this as only suggestive. 
The orbital eccentricity is subject to a positive bias in such fits, because e is positive definite, and experience 
has shown that indications of a small nonzero eccentricity usually disappear after more velocity data are 
obtained. 



- 5 - 



1.005 
J 1.000 ■ 
? 0.995 
| 0.990 
* 0.985 



UT 2006-09-01, Nickel 1m Z, o, = 0.0014 




0.80 



0.85 



0.90 0.95 
HJD - 2453979 



1.00 



1.05 



UT 2006-09-10, Nickel 1m Z, o, = 0.0017 




0.75 



0.80 



0.85 0.90 
HJD - 2453988 



0.95 



1.00 



UT 2006-09-18, FLWO 1.2m z, o, = 0.0016 





0.65 



0.70 



0.75 0.80 

HJD - 2453997 



0.85 



0.90 



UT 2006-09-18, Nickel 1m Z, o, = 0.0019 
*• . • 





0.65 



0.70 



0.75 



0.80 
HJD - 2453997 



0.85 



0.90 



Fig. 1. — Relative photometry of HAT-P-1. The residuals (observed— calculated) are plotted 
beneath the data. 
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Fig. 2. — Relative photometry of HAT-P-1. The residuals (observed— calculated) are plotted 
beneath the data. The bottom panel is a composite light curve created from all of the z and 
Z band data, by subtracting the best-fitting value of T c from the time stamps of each light 
curve and then averaging into 1 min bins. 
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The most natural parameters one would like to know are R+ and R p , but these pa- 
rameters cannot be determined independently from the data. The relevant parameters that 
can be determined are R p /R* and R*/a, where a is the orbital semimajor axis. The only 
property intrinsic to the star that follows directly from the photometric data is the mean 
stellar density (see, e.g., Seager & Mallen-Ornelas 2003): 



In order to determine R* and R p independently, one must have an external estimate of R+, 
or M*, or some combination of R* and M* besides p*. This external estimate may come 
from supplementary observations such as the stellar angular diameter and parallax (see, 
e.g., Baines et al. 2007), or from the interpretation of the stellar spectrum with theoretical 
model atmospheres and evolutionary tracks. The comparison with theory can be facilitated 
by the estimate of p* that comes from the photometry, because the photometric estimate is 
often more precise than the traditional spectroscopic gravity indicator, log g (see, e.g., Pont 
et al. 2007, Sozzetti et al. 2007, Holman et al. 2007). 

Our approach is to fix M* at a fiducial value, and then determine R* and R p from the 
light curve. The scaling relations R± oc M^ 3 and R p oc M^ 3 may then be used to estimate 
the resulting systematic error due to the uncertainty in the stellar mass. This also makes 
it easy to update the determinations of R+ and R p as our understanding of the host star is 
revised through further observations and analyses. In this case, we assumed M* =1.12 M , 
based on the analysis by Bakos et al. (2007) in which the spectroscopic properties of both 
members of the stellar binary were fitted simultaneously to theoretical isochrones. The 
uncertainty in quoted by Bakos et al. (2007) is 8%, corresponding to a systematic error 
of 2.7% in our determinations of R+ and R p . The planetary mass M p hardly affects the 
photometric model at all, since M p <C M*, but for completeness we used the previously 
determined value M p = 0.53 Mj up . 

To calculate the relative flux as a function of the projected separation of the planet 
and the star, we employed the analytic formulas of Mandel & Agol (2002) to compute the 
integral of the intensity over the unobscured portion of the stellar disk. We assumed the 
limb darkening law to be quadratic, 




(1) 



h 



l-a(l-p)-6(l-p) 2 , 



(2) 



where I is the intensity, and p is the cosine of the angle between the line of sight and 
the normal to the stellar surface. We did not use the "small-planet" approximation. We 
fixed the limb-darkening coefficients at the values calculated and tabulated by Claret (2004) 
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for observations of a star with the observed spectral propertiesjj We also investigated the 
effects of changing the limb-darkening law and fitting for the limb darkening parameters, as 
discussed below. In addition, as mentioned in the previous section, we fitted for the zero 
point and slope of the correlation between the measured flux and the external variable that 
showed the strongest correlation. For the FLWO data, this variable was airmass; for Lick it 
was the column number; and for the Wise data it was the FWHM of stellar images. 

The fitting statistic was 



where /j(obs) is the flux observed at time j and Oj controls the weights of the data points, 
and /j(calc) is the calculated flux. As noted in the previous section, the calculated flux was 
the idealized flux of a transit light curve after subtracting a linear function of a specified 
external variable. 

For the data weights (jj, many investigators use the calculated Poisson noise, or the 
observed standard deviation of the out-of-transit data. Experience has shown that these 
procedures usually result in underestimated uncertainties in the model parameters, because 
they neglect time-correlated errors ("red noise"; see, e.g., Gillon et al. 2006), which are 
almost always significant for ground-based data. In order to derive realistic uncertainties on 
the parameters, it is important for <jj to take red noise into account, at least approximately. 

We did this as follows. The most relevant time scale is ~20 min, the ingress or egress 
duration. First we calculated <Ti, the standard deviation of the unbinned out-of-transit data. 
(The results for each light curve are given in Figs. 1 and 2.) Next we averaged the out-of- 
transit data into 20 min bins, with each bin consisting of N data points, depending on the 
cadence. Then we calculated the standard deviation, o^. In the absence of red noise, we 
would observe <Jn = ffi/v~N, but in practice <jat is larger than a\/\fN by some factor f3. 
Therefore, we set the data weights equal to /3 a\. The exact choice of averaging time did not 
matter much. In the end, we used the mean value of (3 over averaging times ranging from 15 
to 25 minutes. Typically we found (3 ~ 2, depending on the telescope and sky conditions^ 



3 Specifically, we used the tables for ATLAS models, interpolating for T e g = 5975 K, logg = 4.45 (cgs), 
log [M/H]= 0.1 and Vt — 2.0 km s _1 . For the z band, a — 0.18 and b = 0.34. We also used these values for 
the Z band, finding it to provide a good fit. For the / band, a — 0.22 and b — 0.34. 

4 This procedure effectively increases the error bar of each measurement and results in a minimum value 
of X 2 /Ndof that is smaller than unity. It is equivalent to setting crj at the value that produce xV-^dof = 1 
but then using Ax 2 = f3 2 instead of Ax 2 = 1 to define the 68% confidence limit. 




(3) 
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In all cases, to solve for the a posteriori probability distributions of the model param- 
eters, we used a Markov Chain Monte Carlo algorithm [see, e.g., Tegmark et al. (2004) 
for applications to cosmological data, Ford (2005) for radial-velocity data, and Holman et 
al. (2006) or Burke et al. (2007) for a similar approach to transit fitting]. We ensured that 
the Gelman & Rubin (1992) R statistic was within 0.5% of unity, a sign of good mixing and 
convergence. For each parameter, we took the median value of the distribution to be our 
best estimate, and the standard deviation as the 1 a uncertainty. 

4. Results 

The results are given in Tables 2 and 3. The first of these tables gives the planetary, 

1 /3 1 /3 

stellar, and orbital parameters, with the fundamental parameters R+/MJ , R p /M p , and i 
listed first. For the parameters that depend on the choice of M* (namely, R*, R p , a, and p p ), 
we have accounted for the systematic error due to the 8% uncertainty in M*. Table 3 gives 
the measured transit times. 



4.1. Planetary, Stellar, and Orbital Parameters 

We find the stellar radius to be R± = 1.115 ± 0.043 R@, and the planetary radius 
to be R p = 1.203 ± 0.051 -Rj up . The statistical error is comparable to the systematic error 
resulting from the covariance with the stellar mass, implying that there is still some scope for 
improvement through additional high-precision photometry. Our value for the stellar radius 
agrees well with the value R± = 1.15^o;o° R & determined by Bakos et al. (2007). Those 
authors estimated -R* from an analysis of the stellar spectrum — its effective temperature, 
surface gravity, and metallicity — whereas we estimated R* (actually R*/M^ 3 ) by fitting 
the transit light curves. The agreement between these different methods of determining 
the stellar radius is an important consistency check on both analyses. Our value for the 
planetary radius is 1.5 a' smaller than the previously determined value R p = 1.36^09 -Rj up , 
where a' is the quadrature sum of the statistical errors of the two estimates. Thus, we have 
revised the planetary radius downward and we have improved the measurement precision by 
a factor of 2. 

We performed a number of additional optimizations to check on the sensitivity of the 
results to the choice of limb darkening function. We tried replacing the quadratic law with 
either a linear law or the four-parameter "nonlinear" law of Claret (2000). For the quadratic 
law, we tried replacing the ATLAS-based coefficients with the PHOENIX-based coefficients 
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of Claret (2004). In none of these cases did the optimized value of R p change by more than 
0.5% relative to the value presented in Table 2. For the case of the linear law, we tried fitting 
for the limb darkening coefficient rather than fixing it at the prescribed value. In that case, 
R p increased by 1.8%, which is still small in comparison to the quoted error. (We found 
that the present data are unable to constrain meaningfully more than one limb-darkening 
parameter.) We conclude that the systematic error due to the choice of limb darkening law 
is probably ~1%. 

Also given in Table 2 are some results reproduced from Bakos et al. (2007) for conve- 
nience, as well as some useful derived quantities. Among these quantities are the impact 
parameter b, defined as a cos i/R* (where a is the semimajor axis), the radius ratio R p /R+, 
the fractions a/R* and a/R p , and the stellar mean density p+, which (as mentioned above) 
do not depend on our choice for M*. We used the previous measurement of the velocity 
semiamplitude of the spectroscopic orbit, K± = 60.3 ± 2.1 m s -1 , to calculate the planetary 
surface gravity, which is also independent of the undetermined stellar properties (see, e.g., 
Southworth et al. 2007, Winn et al. 2007a). The results for a and the planetary mean density 
p p do depend on the choice of stellar mass, and in those cases the quoted errors have been 
enlarged appropriately to take this extra source of uncertainty into account. For convenience 
in planning future observations, we give the calculated values of the full transit duration (the 
time between first and fourth contact, try — tj), and the partial transit duration (the time 
between first and second contact, or between third and fourth contact) E 




Table 3 gives the transit times measured from our data. We have used these times to 
calculate a photometric ephemeris for this system. Using only our 7 new measurements, we 
fitted a linear function of transit epoch E, 



The fit had xV-^dof = 1.6 and iVdof = 5, suggesting that either the period is not exactly 
constant, or that the transit time uncertainties have been underestimated. Because one 
would prefer to have an ephemeris with conservative error estimates for planning future 
observations, we rescaled the measurement errors by y/1.6 and re-fitted the ephemeris, finding 
T c (0) = 2453997.79258(29) [HJD] and P = 4.46543(14) days. The numbers in parentheses 




4.2. 



Transit Times 



T C (E) = T c (0) + EP 



(4) 



5 Although the partial transit duration is listed as t\\ — t\ in Table 1, all of the results in Table 1 are based 
on the entire light curves, including both ingress and egress data. Our model assumes tjj —t\ = try — tm . 
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indicate the 1 a uncertainty in the final two digits. Our derived period agrees with the 
value 4.465290(90) days determined by Bakos et al. (2007), within the respective 1 a limits. 
Figure 3 is the O— C (observed minus calculated) diagram for the transit times. 

5. Summary and Discussion 

We have presented new photometry of HAT-P-1 spanning the times of transit of its 
close-in giant planet. The photometry improves the precision with which the stellar and 
planetary radii are known by a factor of 2, and places the measurements on a more robust 
footing by determining the stellar mean density directly from the photometric data. We 
have also updated the transit ephemeris, to help in planning future observations. 

Although the revised planetary radius is somewhat lower than the previously determined 
value, the planet HAT-P-lb is still among the largest and least dense of the known transiting 
exoplanets. Its mean density (0.376 ± 0.031 g cm -3 ) is comparable to that of the famously 
oversized planet HD 209458b (0.35 g cm~ 3 ; Knutson et al. 2007b). A third planet with 
a comparably low mean density is WASP-lb (0.36 g cm -3 ; Collier Cameron et al. 2007; 
Charbonneau et al. 2007; Shporer et al. 2007). 

There is an extensive literature on the interpretation of exoplanetary radii, and in partic- 
ular on the subject of these apparently "bloated" planets. This term refers to the apparent 
conflict (of order 10-20%) between the measured planetary radii and the calculated radii 
using simple structural models for hydrogen/helium planets of the appropriate mass, tem- 
perature, age, and degree of external heating by the parent star. Many mechanisms have 
been proposed to sustain hotter gaseous envelopes and therefore larger planets: the efficient 
delivery of heat from the star to the planetary interior (Guillot & Showman 2002; Show- 
man & Guillot 2002); the production of internal heat by tidal interactions (Bodenheimer et 
al. 2003; Winn & Holman 2005); and the trapping of internal heat by enhanced atmospheric 
opacities (Burrows et al. 2007) or inhibited convection (Chabrier & Baraffe 2007). Only the 
tidal mechanisms have been specific or predictive enough to be ruled out; the obliquity-tide 
theory of Winn & Holman (2005) did not withstand more detailed dynamical analysis (Lev- 
rard et al. 2007, Fabrycky et al. 2007), and the eccentricity-tide mechanism of Bodenheimer 
et al. (2003) does not seem to be operative in the case for which it was invented, HD 209458b 
(Laughlin et al. 2005). Which (if any) of the other mechanisms contribute to the observed 
radii of transiting exoplanets is not yet clear. 

However, for HAT-P-lb, this issue may be a red herring. We find that there is no 
strong conflict with structural models at this point, as long as the planet does not have 
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Fig. 3. — Transit timing residuals for HAT-P-lb. The calculated times, using the ephemeris 
derived in § 4.2, have been subtracted from the observed times. The filled symbols represent 
observations of complete transits. Open symbols represent observations of partial transits. 
The diamond represents the previous observation by Bakos et al. (2007), which was not used 
in the fit but which agrees well with the computed ephemeris. 
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a massive core of heavy elements. Burrows et al. (2007) recently computed models for 
many of the transiting exoplanets including HAT-P-lb in particular, taking into account the 
appropriate planetary mass, orbital distance, stellar luminosity, stellar spectrum, and stellar 
age. Assuming a planet of solar composition, an atmosphere of standard solar-composition 
opacity, and no dense heavy-element core, they calculated R p = 1.18-1.22 i?j up over the 
plausible age range 3.5 ±1.0 Gyr. This range of calculated values for R p comfortably overlaps 
the 1 a range in our measured value, 1.203 ± 0.051 Rj up . 

Fortney, Marley, & Barnes (2007) have also provided theoretical estimates of exoplane- 
tary radii over a wide range of conditions, although they are not specifically tuned for any 
particular cases of the known exoplanets. For a coreless H/He planet with mass 0.5 M Jup 
orbiting a 3-Gyr-old solar-luminosity star at a distance of 0.045 AU (and thereby receiving 
nearly the same flux as HAT-P-1, which orbits a ~1.5 L & star at a distance of 0.055 AU), 
Fortney et al. (2007) predict a planetary radius R p = 1.12 i?j up . This is smaller than the 
value computed by Burrows et al. (2007), and at least part of the reason for the difference 
is that Fortney et al. (2007) did not account for the "transit radius effect": the enlarged 
size of the transit-measured radius relative to the r = 2/3 photosphere that is usually taken 
to be the "radius" by theoreticians. This effect amounts to a few per cent in the planetary 
radius (see also Burrows et al. 2003). Assuming that this effect adds between 0% and 5% to 
the calculated radius, the difference between the calculated radius of R p and our measured 
value is 0.5-1.6 a, i.e., not very significant. 

We conclude that the present data are consistent with current models of coreless, solar- 
composition, strongly irradiated giant planets. Bakos et al. (2007) estimated a stellar metal- 
licity of Z = 0.025, i.e., comparable to the Sun, and hence the inference of a small or absent 
core is broadly consistent with the core-metallicity relation proposed by Guillot et al. (2006). 
Of course, there are many other possibilities that are also consistent with the data, such as 
a planet with a dense core that also has either an extra source of atmospheric opacity or 
an extra source of internal heat. As of now there is no way to distinguish among these 
possibilities. 

As discussed by Burrows et al. (2007), it is becoming clear that there are many de- 
terminants of planetary radii. By considering the entire ensemble of exoplanets one can 
fully appreciate the strengths and weaknesses of theoretical models, and possibly obtain 
clues about interesting processes that may have been overlooked. This requires not only the 
discovery of new transiting systems, but also high-precision determinations of the system 
parameters, such as the present study. 

We thank Debra Fischer and Geoff Marcy for helpful discussions, and John South- 
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Table 1. Photometry of HAT-P-1 



Telescope 


Filter 


Heliocentric Julian Date 


Relative flux 


FLWO 


z 


2453997.69528 


0.9988 


FLWO 


z 


2453997.69560 


1.0006 


FLWO 


z 


2453997.69591 


1.0031 


FLWO 


z 


2453997.69621 


1.0009 


Lick 


Z 


2453997.91750 


1.0013 


Lick 


Z 


2453997.91781 


0.9965 


Lick 


Z 


2453997.91811 


1.0016 


Wise 


I 


2454069.33838 


1.0014 


Wise 


I 


2454069.33894 


0.9970 


Wise 


I 


2454069.33950 


1.0000 



Note. - The time stamps represent the Heliocentric Julian 
Date at the time of mid-exposure. We intend for this Table to 
appear in entirety in the electronic version of the journal. Ex- 
cerpts are shown here to illustrate its format. The data are also 
available from the authors upon request. 
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Table 2. System Parameters of HAT-P-lb 



Parameter 


A 7" 1 

Value 


Uncertainty 


(R lRrM(M 11 12 Mn) 1/3 


1.115 


0.034 


fp /p. \I(M /I 19 M^W 3 


1 903 


U.U'iO 


7 [ripen 


86.22 


0.24 




1 19 


n dq 


A/f [Mt 1 




n 04 

U.U1 


A/jdlo r*l "l"\7" Qf^TYl 1 TYlTll 11 f\ & TV \\W Q "M 
V ClVj KjL by QCllllCllllUll \j « ± V ^ ill O I 


fin s 


9 1 


1 IvKlT O T~\ DTI r\ *H /— ' H OTTO 

vjiuiiai peiioci, r [ciaysj 


4 4fi^9Q 


n nnnno 
u.uuuuy 


D [D 1 


1. 110 




-R» -Rjup 

if L 1 J 


1.203 


0.051 


-Rp/ -R* 


0.11094 


0.00082 


Semimajor axis, a [AU] 


0.0551 


0.0015 


6 = a cos i/R± 


0.701 


0.023 


a/R* 


10.64 


0.32 


a/ R p 


95.9 


3.5 


t w - ti [hr] 


2.779 


0.032 


in - h [hr] 


0.508 


0.035 


p* [g cm -3 ] 


1.14 


0.10 


p p [g cm" 3 ] 


0.376 


0.031 


GM P / Rp [cm s- 2 ] 


904.5 


66.1 



Note. — This table has three sections, divided by horizontal 
lines. The top section lists the parameters that were estimated 
by fitting the new photometric data, as explained in § 3. The 
orbital eccentricity was assumed to be exactly zero. The mid- 
dle section lists some parameters from Bakos et al. (2007), re- 
produced here for convenience. The bottom section lists some 
interesting parameters that can be derived from the parameters 
in the first two sections. 
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Table 3. Mid-transit times of HAT-P-1 



Observatory 


Epoch 


Mid-transit time 


Uncertainty 




E 


[HJD] 


[days] 


Lick 


-4 


2453979.92848 


0.00069 


Lick 


-2 


2453988.86197 


0.00076 


FLWO 





2453997.79200 


0.00054 


Lick 





2453997.79348 


0.00047 


FLWO 


2 


2454006.72326 


0.00059 


FLWO 


3 


2454011.18837 


0.00107 


Wise 


16 


2454069.23795 


0.00290 



Note. - - Based on these new measurements, we derived 
a transit ephemeris T C (E) = T c (0) + EP with T c (0) = 
2453997.79258(29) [HJD] and P = 4.46543(14) days, where 
the numbers in parentheses indicate the la uncertainty in 
the final two digits. We note that Bakos et al. (2007) de- 
rived a more precise period based on observations over 217 
cycles, P = 4.465290(90) days. 



